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Abstract. The BCS and HFB theories which can accommodate the pairing correlations in 
the ground states of atomic nuclei are presented. As an application of the pairing theories , 
we investigate the spatial extension of weakly bound Ne and C isotopes by taking into account 
the pairing correlation with the Hartree-Fock-Bogoliubov (HFB) method and a 3-body model, 
respectively. We show that the odd-even staggering in the reaction cross sections of 30 > 31 > 32 Ne 
and 14 . 15 . 16 C are successfully reproduced, and thus the staggering can be attributed to the unique 
role of pairing correlations in nuclei far from the stability line. A correlation between a one- 
neutron separation energy and the anti-halo effect is demonstrated for s- and p-waves using the 
HFB wave functions. We also propose effective density-dependent pairing interactions which 
reproduce both the neutron-neutron (nn) scattering length at zero density and the neutron 
pairing gap in uniform matter. Then, we apply these interactions to study pairing gaps in 
semi-magic finite nuclei, such as Ca, Ni, Sn and Pb isotopic chains. 



1. Introduction 

It has been known that the pairing correlations play an important role in finite and also infinite 
nuclear systems. Just after the BCS theory was proposed as a fundamental theory of metallic 
superconductor, Bohr, Mottelson and Pines proposed a possible analogy of superfluidity in 
atomic nuclei in 1958 [1]. The most prominent evidence for pairing correlation in nuclei is 
found in the odd-even staggering in binding energies and the energy gap in the excitation 
spectrum of even-even nuclei in contrast to a compressed quasi-particle spectrum in odd-A nuclei 
[1, 2, 3]. There are also dynamical effects of pairing correlations seen in the moment of inertia 
associated with nuclear rotation and large amplitude collective motion [3, 4, 5]. The Hartree- 
Fock (HF)+BCS method and Hartree-Fock-Bogoliubov (HFB) method have been commonly 
used to study the ground state properties of superfluid nuclei in a broad mass region [6, 7, 8, 9]. 

As new phenomena in nuclei near the neutron-drip line, large odd-even staggering (OES) 
phenomena have been revealed experimentally in reaction cross sections of the isotopes 14 > 15 > 16 C 

[10], 18 > 19 ' 20 C [11], 28 < 29 < 30 Ne[12], 30^1,32 Ne [ 12 ] > and 36,37,38 Mg In Ref wg haye argued 

that the pairing correlations play an essential role in these OES. That is, the OES in reaction 
cross sections is intimately related to the so called pairing anti-halo effect discussed in Ref. [15]. 
In this lecture , we discuss the pairing correlations close to the zero energy by a Hartree-Fock 
Bogoliubov (HFB) method. This problem is also related with the superfluidity of neutron gases 
in the outer crust of neutron stars [16]. 



Recently, we proposed new types of density dependent contact pairing interaction, which 
reproduce pairing gaps in a wide range of nuclear mass table [9]. We discussed also the relation 
between the proposed paring interactions and the pairing gaps in symmetric and neutron matters 
obtained by a microscopic treatment based on the nucleon-nucleon interaction. We will show 
the necessity of the isovector type pairing interaction on top of the isoscalar term to reproduce 
systematically nuclear empirical pairing gaps. 

This lecture note is organized as follows. Section 2 gives basic formulas of BCS and HFB 
theories. In section 3, we will discuss the odd-even staggering in the reaction cross sections in 
nuclei near the neutron drip line. Section 4 is devoted to study the isospin dependent pairing 
interaction for the study of pairing correlations in semi-magic nuclei. 



2. BCS and HFB theories 

A very important generalization of the HF theory is to accommodate pairing interaction on 
top of the mean field. The general theoretical framework with the pairing correlations for the 
single-particle orbitals is called Hartree-Fock-Bogoliubov (HFB) equations, that is also called 
the Bogoliubov-deGenne equations in condensed matter physics. A simpler version of HFB is 
called the BCS theory which has been often employed in nuclear physics. In dealing with even- 
even nuclei, the HF equations are invariant under time reversal. This implies that each orbital i 
has its time-reversed partner i, and the two orbits have the same single-particle energy £j. The 
basic BCS Ansatz for the many-particle wave function is 

|BCS) = n(« i + Ui a{o|)|0) 

i>o 

where aj is the particle creation operator acting on the HF vacuum |0). The parameters U{ and Vi 
will be determined by minimizing the expectation value of the Hamiltonian. The normalization 
of the state requires 

\ui\ 2 + \vi\ 2 = 1 
The BCS state (1) is further rewritten to be 

|BCS}ocexp(^^aJa})|0). (3) 

The BCS state is not an eigenstate of the particle number. In condensed matter physics, this 
does not cause any serious problem since the number of particles is close to the Avogadro number 
~10 23 . On the other hand, in nuclear physics, the number of nucleons is in the order of at most 
200 so that we have to take care of the number conservation in some way. One possible way is 
to introduce a Lagrange multiplier term in the Hamiltonian 

H' = H — XN. (4) 

Then the particle number expectation value can be fixed to the desired value N 

(BCS|iV|BCS) = 2 ^2 v f = N (5) 

i>0 

on average. The quasi-particles are introduced by a unitary transformation, so called the 
Bogoliubov-Valatin transformation 

a J = UiCb\ — Vicq 

cq = Uidj + Via} (6) 



(1) 



(2) 



where Oj is the physical (or bare) annihilation operator. For the unitary transformation, the 
quasi-particles preserve the anti-commutation relations 

{a.i,ai)} = Si tj , {a i: aj} = {a\,at} = 0. (7) 

The BCS vacuum is defined by an equation 

ai|BCS) = 0. (8) 

We can prove that the definition (8) is equivalent to the BCS vacuum (1). Firstly, the definition 
(8) can read 

ibcs)= n «iio> = n «i«ii°>- (9) 

iGall i>0 

By performing the transformation (6) from the quasi-particles to the bare particles in Eq. (9), 
we obtain 

|BCS) = ]J{u 2 ajcq + Ujt>;(l - a|a,) - u^ata^ + vfa\at}\0) 

i>0 

= II Vi{ui + v ia \a^}\0) a Y[{ Ul + ^aja|}|0). (10) 

i>0 i>0 

The hamiltonian density can be evaluated for the BCS state as 

<BCS|iT|BCS) = £( £i - \)v 2 i - J2 ViijjUiViUjVj (11) 

i i,j>0 

The variation of the Hamiltonian density with respect to the Vi (or equivalently Ui ) gives a BCS 
gap equation to be solved; 

A, = E 2jl P 1V iifj u ^ (12) 

j>0 



with 



u, = - 1 + 



^/( £i -A)2 + A? 



"f = f 1- / (13) 

The quasi-particle energy Ei is evaluated by a commutator relation 

= (14) 

to be 

Ei = \J (ei - A) 2 + A 2 . (15) 

The BCS theory is generalized to HFB theory by using the following unitary transformation. 
The Hartree-Fock-Bogoliubov (HFB) transformation for quasi-particle can read 

a l X "v"J '•<./"., 

j 



where Uy and vij are the upper and the lower components of BCS transformation, respectively. 
The normalization condition for the Uij and factors is given by 

£u?-+4 = l. (17) 

3 

The HFB model can be generalized to the coordinate space representation as [17] 

a\ = J dr{u t (r)ip\r) + v*(r)iP(r)} 

a t = J dr{<(r)^(r) + ^(r)^(r)} (18) 

where ip^(r) and ip(r) are the creation and the annihilation operators of bare particle in the 
coordinate space representation and obey the anti-commutator relation 

{^(r),^ f (r')} = 5(r-r'). (19) 
The orthonormal conditions of Ui (r) and Vi (r) functions are given by 



dr{u*(r) Uj (r) + v*(r) Vj (r)} = S id 
^{u*(r) Ui (r') + v*(r) Vi (r')} = <5(r - r'). (20) 

i 

The density and the pair density (abnormal density) are obtained as 

p(r,r') = (HFB|^(r)Y>(r')|HFB) = ^« t (r)<(r') (21) 

i 

k{y, r') = (HFB|^(r)Y>(r')|HFB) = ]T v^r)^) (22) 

i 

where the HFB vacuum is defined by 

ai |HFB) = 0. (23) 
The hamiltonian density for the constrained Hamiltonian H' (4) is expressed to be 

(HFB | H' | HFB) = Tr{(T - X)p} + ^Tr{Vpp} + ^{Vkk} (24) 

where T and V are the kinetic density and the hamiltonian density due to a two-body interaction, 
respectively. 

The HFB equations are obtained by the variation of the hamiltonian density in Eq. (24) with 
respect to the itj(r) and Vi(r) functions to be 

T + Vhf-X A(r) ^ ( Ui (r) \ _ „ ( Ui (r 

A(r) 

where Vhf is the HF potential, 



A(r) -T - V HF + A A Vi(r) J ~ Ei [ v]{v) > " (25) 



Vhf = Tr{V P } (26) 



and A(r) is the pairing gap potential, 



A(r) = Tr{VK}, (27) 

respectively. In Eq. (25), u«(r) and Vi(r) are radial wave functions and could have different 
shapes from the HF single particle wave function, which is the eigenstate of HF hamiltonian 
huF = T — Vhf- In contrast, in the BCS model, the Ui and vi are factors and simply multiplied 
to the HF single particle wave function fi by the Bogoliubov-Valatin transformation (6). 

A question is whether the HFB is necessary for realistic framework to treat the pairing 
correlations in nuclei. When the active orbitals are well bound as in nuclei along the valley of 
stability, the single-particle wave functions are not influenced by the pairing field and the BCS 
model could be entirely adequate. This might not be the case in nuclei near the drip lines, 
where the pairing correlations would make the difference between bound and unbound orbitals. 
In a recent paper, we showed that a major part of the effects of HFB could be included by a 
perturbative treatment of the BCS wave functions even in nuclei far from the valley of stability 
[18]. 



3. Halo and Pairing correlations 

The pairing correlations play a crucial role to develop the halo structure of loosely bound nuclei. 
At the same time, the pairing will act to prevent the divergence of very loosely bound nuclcons 
in the mean field. Let us discuss how the single-particle wave function will expand its tail in 
the mean field potential and how the pairing correlations will prevent the exponential growth of 
expansion to be infinity. In the mean field model, the single particle wave function 

^m(r) = Vim(r)Y lm (r) = ^-Y lm (r) (28) 
is calculated with a Schrodinger equation 

/ h 2 d 2 l(l + i)h 2 , . \ 

(-^ + -4^ + ^< r »- j ')" ,<r)=a (29) 

The asymptotic behavior in the limit r — >■ oo with Vhf —> is given by the modified spherical 
Bessel function and behaves 

ui=o(r) — >■ exp(— nr) in the limit r — > oo (30) 

for the angular momentum I = case. Then the mean square radius for s— wave is evaluated to 
be 

, 2 v _ = J r 2 ui =0 (r) 2 dr ^ l_ = % 2 _^ ^ 
_0 / ui =0 (r) 2 dr k 2 2m\e\ 

which diverges in the limit of very loosely bound case |e| — > 0. It is also shown in the case of 
p— wave case, 1 = 1, that the mean square radius also diverges but slightly slowly proportional 
to l/-y/je[ [19]. On the other hand, the wave function with I > 2 does not show any divergence 
of mean square radius in the limit of |e| — > [20]. 

In contrast, the upper component of a HFB wave function, which is relevant to the density 
distribution, behaves [21] 

v(r) oc exp(-/3r), (32) 
where j3 is proportional to the square root of the quasi-particle energy E, 



/ 3 = \Tt(^-A), (33) 



A being the chemical potential. If we evaluate the quasi-particle energy in the BCS 
approximation or HFB with canonical basis, it is given as 



£ = v /( £ -A)2 + A2, (34) 

where A is the pairing gap. For a weakly bound single-particle state with e ~ and A ~ 0, the 
asymptotic behavior of the wave function v(r) is therefore determined by the gap parameter as, 



v(r) oc exp | — y -^-A • r ] . (35) 

The radius of the HFB wave function will then be given in the limit of small separation energy 
|e| — >■ as 

, 2^ Jr 2 \v(r)\ 2 dr 1 fi 2 

{r >HFB = j\ v{rWdr <* jga ^A" (36) 

As we show, the gap parameter A stays finite even in the zero energy limit of e with a density 
dependent pairing interaction. Thus the extremely large extension of a halo wave function in 
the HF field will be reduced substantially by the pairing field and the root-mean-square (rms) 
radius of the HFB wave function will not diverge. This is called the anti-halo effect due to the 
pairing correlations [15]. 

We now numerically carry out mean-field calculations with a Woods-Saxon (WS) potential 
and also HFB calculations using the single-particle wave functions in the WS potential. Most of 
following materials of this section are taken from refs. [14]. As examples of s-wave and p-wave 
states, we choose the 2s 1 / 2 state in 23 and 2p 3 / 2 state in 31 Ne, respectively. Although 31 Ne 
is most likely a deformed nucleus[22, 23], for simplicity we assume a spherical Woods-Saxon 
mean-field potential. Notice that a Woods-Saxon potential with a large diffuseness parameter a 
yields the 2p 3 / 2 state which is lower in energy than the I/7/2 state, as was shown in Ref. [24]. 
We use a similar potential with a=0.75 fm as in Ref. [24] for 31 Ne, while that in Ref. [25] 
for 23 O. For the HFB calculations, we use the density-dependent contact pairing interaction of 
surface type, in which the parameters are adjusted in order to reproduce the empirical neutron 
pairing gap for 30 Ne [26]. While we fix the Woods-Saxon potential for the mean-field part, the 
pairing potential is obtained self-consistently with the contact interaction. 

The left top panel of Fig. 1 shows the mean square radius of the 2si/ 2 state for 23 0, while 
right top panel of Fig. 1 shows the mean square radius of the 2p 3 / 2 state for 31 Ne. In order to 
investigate the dependence on the single-particle energy, we vary the depth of the Woods-Saxon 
wells for the Si/ 2 and j> 3 / 2 states for 23 and 31 Ne, respectively. The dashed lines are obtained 
with the single-particle wave functions, while the solid lines are obtained with the wave function 
for the canonical basis in the HFB calculations. One can see an extremely large increase of the 
radius of the WS wave function for both the s-wave and p-wave states in the limit of e\ys 0. 
In contrast, the HFB wave functions show only a small increase of radius even in the limit of 
e ws — * 0. This feature remains the same even when the contribution of the other orbits are 
taken into account, as shown in the middle panel of Figs. 1. Due to the pairing effect in the 
continuum, the HFB calculations yield a larger radius than the HF calculations for the cases 
of ews < — 1 MeV. On the other hand, in the case of —1 MeV < e\ys < MeV the HF wave 
function (equivalently one quasi-particle wave function in HFB) extends largely, while the HFB 
wave function does not get much extension due to the pairing anti-halo effect. 

In the bottom panel of Figs. 1, the average pairing gaps are shown as a function of the single 
particle energy ews- It is seen that the average pairing gap increases as ews approaches zero. 
This is due to the fact that the paring field couples with the extended wave functions of weakly 
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Figure 1. (Color online) (left panel) The mean square radii and the average paring gap as 
a function of the single particle energy ews in a Woods-Saxon mean-field potential. The top 
panel shows the mean square radius of the 2si/2 wave function with and without the pairing 
correlation, denoted by HFB and Woods-Saxon, respectively. The middle panel shows the rms 
radii for 22 (the dotted line), 23 (the dashed line), and 24 (the solid line), obtained with the 
Hartree-Fock ( 22 and 23 0) and the Hartree-Fock-Bogoliubov ( 24 0) calculations. The bottom 
panel shows the results of the HFB calculations for the average paring gap of 24 O. (right panel) 
Same as the left panel, but for the 2p 3/ / 2 state and for 30 > 31 > 32 Ne isotopes. 



bound nucleons in the self-consistent calculations. That is, the pairing field is extended as the 
wave functions do and becomes larger for a loosely bound system. 

Let us now calculate the reaction cross sections for the 30 > 31 > 32 Ne isotopes and discuss the 
role of pairing anti-halo effect. To this end, we use the Glauber theory, in which we adopt the 
prescription in Refs. [24, 27] in order to take into account the effect beyond the optical limit 
approximation. The left panel of Fig. 2 shows the reaction cross sections of the 30 > 31 > 32 Ne nuclei 
on a 12 C target at 240 MeV/nucleon. We use the target density given in Ref. [28] and the profile 
function for the nucleon-nucleon scattering given in Ref. [29]. In order to evaluate the phase 
shift function, we use the two-dimensional Fourier transform technique [30]. The cross sections 
an shown in Fig. 2 are calculated by using projectile densities constrained to two different 
separation energies of the 2p 3/ / 2 neutron state. The dashed line with triangles is obtained using 
the wave functions with the separation energy S n ( 31 Ne) =|ews|=l-49 MeV, while the solid line 
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Figure 2. (Color online) (left panel) Reaction cross sections of Ne isotopes on a 12 C target at 
-£'iab = 240 MeV/A. The cross sections are calculated with the Glauber theory with HF and HFB 
densities. The solid line with the filled squares shows the results of 5 n ( 31 Ne) =|e(2p 3/2 ) |=0.32 
MeV, while the dashed line with the open triangles is obtained for S' n ( 31 Ne) =|e(2p 3 / 2 )|=l-49 
MeV. The experimental data are taken from Ref. [12]. (right panel) Reaction cross sections of C 
isotopes on a 12 C target at i?i a b=83 MeV/A. The cross sections are calculated with the Glauber 
theory with three-body model densities. The experimental data are taken from Ref. [10]. 



with squares is calculated with the wave functions of 5' ra ( 31 Ne)=|ews |=0.32 MeV. The empirical 
separation energy of 31 Ne has a large ambiguity with SW=0.29±1.64 MeV[31]. The cross section 
ctr of 30 Ne is already much larger than the systematic values of Ne isotopes with A<30. On 
top of that, we can see a clear odd-even staggering in the results with the smaller separation 
energy, S' n =0.32 MeV, as much as in the experimental data, while almost no staggering is seen 
in the case of the larger separation energy, S n =lA9 MeV. This difference is easily understood by 
looking at the anti-halo effect for |e\ys| < lMeV shown in the middle panel of Fig. 1. Recently, 
the effect of deformation of neutron-rich Ne isotopes on reaction cross sections was evaluated 
using a deformed Woods-Saxon model [32, 33]. It was shown that the deformation is large as 
much as /?2 ~ 0.42 in 31 Ne and enhances the reaction cross section by about 5%. However, the 
calculated results did not show any significant odd-even staggering in between 28 Ne and 32 Ne 
[32]. 

We investigate next neutron-rich C isotopes. We particularly study the 16 C nucleus using 
a three-body model given in Ref. [34]. In this case, the valence neutron in 15 C occupies the 
2s 1 /2 level at ews = —1-21 MeV, while 16 C is an admixture of mainly the ^s^) 2 an d (ld5/2) 2 
configurations. Assuming the set D given in Ref. [34] for the parameters of the Woods-Saxon 
and the density distribution for 14 C given in Ref. [35], the rms radii are estimated to be 2.53, 
2.90, and 2.81 fm for 14 C, 15 C, and 16 C, respectively. The corresponding reaction cross sections 
<jr calculated with the Glauber theory are shown in the right panel of Fig. 2. The calculation 
well reproduces the experimental odd-even staggering of the reaction cross sections, that is a 
clear manifestation of the pairing anti-halo effect. 

In order to quantify the OES of reaction cross sections, we introduce the staggering parameter 
defined by 

, N W-4 + 1) - 2a R (A) + a R (A - 1) 

73 = (-) o ' ( 67 ) 



where ctr(A) is the reaction cross section of a nucleus with mass number A. We can define the 
same quantity also for rms radii. Notice that this staggering parameter is similar to the one 
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Figure 3. (Color online) The staggering parameter 73 of reaction cross sections denned by 
Eq. (37) for the Ne isotopes with the 12 C target at £'=240 MeV/nucleon. This is plotted as 
a function of the neutron separation energy S n of the odd-yl nuclei. The experimental data 
for the reaction cross sections are taken from Ref. [12], while the empirical separation energies 
are taken from Refs. [36, 31]. The dashed line is the calculated staggering parameter for the 
30,3l,32]vj e isotopes, assuming that the valence neutron of 31 Ne occupies the 2p3/2 orbit. 



often used for the OES of binding energy, that is, the pairing gap. 

The experimental staggering parameters 73 are plotted in Fig. 3 for Ne isotopes as a function 
of the neutron separation energy for the odd-mass nuclei. We use the experimental reaction 
cross sections given in Ref. [13] while we evaluate the separation energies with the empirical 
binding energies listed in Ref. [36]. For the neutron separation energy for the 31 Ne nucleus, 
we use the value in Ref. [31]. The experimental uncertainties of the staggering parameter are 

obtained as 

, VWgA + 1))* + A{5a R {A)Y + (5a R (A - l)f 

where 5<jr{A) is the experimental uncertainty for the reaction cross section of a nucleus with 
mass number A. The figure also shows by the dashed line the calculated staggering parameter for 
the 30 > 31 > 32 Ne nuclei with the 2p 3 / 2 orbit. One sees that the experimental staggering parameter 
agrees with the calculated value for 30 ' 31 ' 32 Ne nuclei when one assumes that the valence neutron 
in 31 Ne occupies the 2p 3/ / 2 orbit. Furthermore, although the structure of lighter odd-A Ne 
isotopes is not known well, it is interesting to see that the empirical staggering parameters 
closely follow the calculated values for the 2p 3/ / 2 orbit. This may indicate that the low-Z single- 
particle orbits are appreciably mixed in these Ne isotopes due to the deformation effects [22, 23]. 

4. Isospin Dependent Pairing Interaction and Pairing Gaps in Finite Nuclei 

There are mainly two different approaches for a calculation of pairing correlations in finite 
nuclei. The first approach is based on phenomenological pairing interactions whose parameters 
are determined using some selected data [37], while the second approach starts from a bare 



nucleoli- nucleon interaction and eventually includes the effect of phonon coupling [38]. The 
latter approach has shown that the medium polarization reduces the pairing gaps in neutron 
matter while, in symmetric matter, the neutron pairing gaps are much enlarged at low density 
compared to that of the bare calculation. This enhancement takes place especially for neutron 
Fermi momenta /cFn < 0.7 fm _1 . 

We propose effective density-dependent pairing interactions which reproduce both the 
neutron-neutron (nn) scattering length at zero density and the neutron pairing gap in uniform 
matter. In order to simultaneously describe the density dependence of the neutron pairing gap 
for both symmetric and neutron matter, it is necessary to include an isospin dependence in the 
effective pairing interaction [9]. Depending on whether the medium polarization effects on the 
pairing gap given in Ref. [39] are taken into account or not, we invent two different density 
dependent functionals in the pairing interaction. Then, we apply these interactions to study 
pairing gaps in semi-magic finite nuclei, such as Ca, Ni, Sn and Pb isotopic chains. 

The density— dependent pairing interaction can be read as 

Upair(l, 2) = g[p, I] <*(?! " ? 2 ). (39) 

where p is the nuclear density and I is defined as I = (p n — p p ) / p ■ t z [9]. In Ref. [9], an isovector 
dependence in the density-dependent term g is separated to two parts g = g 1 + g 2 . The function 
g 1 is determined to mimic the bare pairing gaps in nuclear matters and the function g 2 takes 
care of the medium polarization effect. The functional form of g 1 is given by 

g 1 ^,/] = i - ui) Vs (p-Y s - ui) Vn (p-Y n , (40) 

\PoJ \PoJ 

where po=0.16 fm~ 3 is the saturation density of symmetric nuclear matter and the functions 
f s (7) and f n (J) are f s (7) = 1 - f n (J) and f n (J) = I = (p n (r) - p P (r))/p(r) ■ r z . The values of 
parameters r] s ,rj Q and powers of density dependence a s ,a n are given elsewhere [9]. 

We perform Hartree-Fock-Bogoliubov (HFB) calculations for semi-magic Calcium, Nickel, 
Tin and Lead isotopes using these density-dependent pairing interactions in Eq. (39) derived 
from a microscopic nucleon-nucleon interaction. Our calculations reproduce well the neutron 
number dependence of experimental data for binding energy, two neutrons separation energy, 
and odd-even mass staggering of these isotopes [40]. Especially the interaction IS+IV Bare 
without the medium polarization effect gives satisfactory results of all the isotopes as is seen 
in Fig. 4. It is clear in the comparison between IS bare and IS+IV bare that the isospin 
dependence of the pairing interaction plays an important role in the pairing gaps of neutron-rich 
nuclei. The isospin dependent pairing interaction is further applied to both even-even and even- 
odd nuclei by using EV8-odd program. The results successfully reproduce the empirical isotope 
and isotone dependent functionals of the odd-even mass differences of several medium-heavy 
and heavy nuclei [41]. Yamagami et al., studied also the quadratic isospin dependent terms of 
pairing interaction [42, 43]. They adopted several Skyrme interactions with different effective 
mass and showed the effective mass dependence of the isospin dependent terms. It is also shown 
in ref. [44] that the effect of Coulomb two-body interaction on the proton pairing gaps can be 
mimicked by the leaner isospin dependent term in the pairing interaction. 

5. Summary 

We have studied the mass radii of Ne isotopes with the Hartree-Fock (HF) and Hartree-Fock- 
Bogoliubov (HFB) methods with a Woods-Saxon potential. The reaction cross sections (Tr were 
calculated using the Glauber theory with these microscopic densities. We have shown that the 
empirical odd-even staggering in the reaction cross sections of neutron-rich Ne isotopes with 



4 

3 

a.fA 



Ca 



HFB IS+IV Bare 

HFB IS+IV Induced 

HFB IS Bare 



Ni 



V 



j^ l ^' l I I V|N - < I l 'f I fj frf 'l 'i I I f' l I I I '\' \ r^H | Hr 



2 

^ c 1 

< 
tO -1 

-2 

4 

3 



2 

< 
tO -1 

-2 



20 



30 



40 



30 



40 



50 



60 



1 1 1 1 1 - 

L Sn J 

Tffiji'rri~f'i 1 1 rVrnfi'i 1 1 1 1 1 1 1 1 1 rjf 


_ | ! ! 1 | 1 1 1 | 1 1 1 | 1 1 . | . ; 

~r Pb ^ 

It i i i | N rfV1'T i i | i i i v |V 


- — "~ ■ — jj, , . . ■ 


- ■ ^* . 




_ *v ' 

1 l l l 1 t t t f r r r [ T T l 1 l 



50 60 70 



80 90 

N 



100 110 120100 



120 



140 

N 



160 



180 



Figure 4. Comparison of the neutron HFB pairing gaps A n with the odd-even mass staggering 
given by the three-point formula . The dotted line shows the results of the pairing 
interactions IS+IV Bare, while the short dashed, and long dashed lines are obtained with the 
pairing interactions IS+IV Screened and IS Bare, respectively. The difference 5(A n ) is defined 
as 5(A n ) = A n (th.) — A n (exp.). All units are given in MeV. 



the mass A=30 ~ 32 is well described by the HFB density and can be considered as a clear 
manifestation of the pairing anti-halo effect associated with a loosely-bound 2p 3 / 2 wave function. 
The index of the odd-even staggering is proposed and applied successfully for the study of the 
reaction cross sections of Ne isotopes to revel the role of the pairing correlations in exotic nuclei. 

A new type of density dependent contact pairing interaction was obtained to reproduce the 
microscopic pairing gaps in symmetric and neutron matter. We performed also HFB calculations 
for semi-magic Calcium, Nickel, Tin and Lead isotopes using these density-dependent pairing 
interactions. Our calculations reproduce well the neutron number dependence of experimental 
data of binding energy, two neutrons separation energy, and odd-even mass staggering of 
these isotopes with the interaction IS+IV Bare. Recently, the isospin dependent interaction 
is extended to introduce the quadratic terms of isospin and applied to a wide region of nuclei 
in the mass table. These result suggests that by introducing the isospin dependent term in the 
pairing interaction, one can construct a global effective pairing interaction which is applicable 
to nuclei in a wide range of the nuclear chart. 
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